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EXPERIMENTAL INVESTIGATION OF THE THERMAL 
BOUNDARY LAYER IN A CONFINED LIQUID 
by Michael J. Kolar 
Lewis Research Center 

SUMMARY 

The boundary layer flow in a confined liquid is studied in the case of constant heat 
flux. Parameters are developed to indicate the presence of both a quasi-steady state and 
boundary layer flow. The conditions under which mathematical similarity may exist in 
the boundary layer temperature and velocity fields are derived, and the coupling between 
the boundary layer and interior flows is discussed. Measured temperature profiles in 
the boundary layer exhibit similarity in a region in which startup and free surface effects 
are negligible. The predicted temperature profiles for a flat plate in an infinite medium 
agree fairly well with the measured profiles in the region where similarity exists. 


INTRODUCTION 

The general problem of a confined liquid subjected to a heat flux has been studied re- 
cently in connection with the heating of cryogenic propellants for rockets. In order to 
design propellant tanks for space vehicles properly, the behavior of the propellant flow 
field should be known. Small scale experiments with noncryogenic liquids have been re- 
ported in references 1 and 2, while larger scale experiments with liquid hydrogen have 
been reported in reference 3. An analytical solution to the momentum and energy equa- 
tions governing a liquid contained in a horizontal cylinder subjected to wall heating is re- 
ported in reference 4. Numerical solutions to the Navier-Stokes equations for liquid 
contained in a vertical parallelepiped with heated walls are given in reference 5, which 
also gives an excellent bibliography on confined heated liquids. These investigations indi- 
cate that, under the heat loads usually experienced by a space vehicle during powered 
flight, a natural- convection boundary layer forms along the propellant tank wall, and an 
axial temperature gradient develops in the fluid bulk. In addition, a layer of saturated 
liquid forms near the surface. The analyses also show that the temperature throughout 



the tank depends very slightly on the radial coordinate. 

When an analytical approach is used to predict the flow field in a rocket propellant 
tank, a number of engineering approximations are usually made to simplify the Navier- 
Stokes equations. The simplifying assumptions usually proceed with the application of 
boundary layer theory taken from steady- state results on vertical flat plates in an infinite 
medium. Therefore, the approximate solutions obtained in propellant heating analyses 
preclude possible transient effects and ignore the coupling that exists between the bound- 
ary layer and the interior, or bulk, flow. Neglecting the transient effects is based on the 
concept that the fluid is in a quasi- steady state, that is, that the flow fields behave as if 
they were governed by the steady-state process at each instant of time. Ignoring the cou- 
pling between the boundary layer and interior flows is generally based on the argument 
that the propellant tank is sufficiently large so that only small changes can occur in the 
interior flow, and therefore, that only small changes can occur in the boundary layer as 
a result of the coupling. In other words, only a weak coupling exists between the bound- 
ary layer and interior flows. Under these conditions, the flow field in a completely con- 
fined fluid is replaced by free- convection flow in an infinite medium. The assumption is 
then made that the temperature distributions in the tank wall and in the propellant are 
such that the conditions of profile similarity for both the thermal and velocity boundary 
layers are satisfied. 

The conditions for the aforementioned assumptions to be valid are neither quantita- 
tively defined nor universally met. Thus, to satisfy the need for a more definitive de- 
scription of completely confined fluids, a study was made to determine criteria under 
which the preceding assumptions are valid. The purpose of this study was twofold: (1) to 
arrive at a set of parameters that define the necessary conditions under which the simpli- 
fying assumptions can be made and (2) to verify experimentally that these conditions can 
be satisfied for a practical situation. 

An analysis was made, which consisted of normalizing the equations which govern 
the flow field in a confined fluid and studying the resulting parameters. Experiments 
were performed in a glass tank, which enclosed a noncryogenic fluid and which was sub- 
jected to a constant wall heat flux. Experimental data were then compared with the ana- 
lytical parameters that were developed. 


GOVERNING PARAMETERS 

The discussion that follows is based on the flow in a two-dimensional rectangular 
cell of liquid subjected to a constant wall heat flux (fig. 1). The fluid is considered to 
have a free surface on the top and thermal insulation at the bottom. It is assumed that 
all fluid properties are constant, and viscous dissipation and work against the body force 
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Figure 1 - Proposed model for enclosed liquid subjected to constant 
wall heat flux. 


are neglected. Density is assumed constant except in the formation of buoyancy terms. 
The coordinate system is taken with the origin at one corner of the cell, x as the ordi- 
nate, and y as the abscissa. The cell is assumed symmetrical about the centerline. 
Only the laminar flow regime is considered. 


Conservation Laws 

The equations expressing conservation of mass, x- and y-components of momentum, 
and energy in an incompressible, two-dimensional flow field may be written as (ref. 6) 

o 

ax ay 

^ + u^ + v^ = l-l 9 p + l /V 2 u 

3r 3x 3y p p 3 x 

3v dv dv 1 3p _2 

3r 3x 3y p 9y 
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field and wall conditions proceeds as though a steady state process exists at each instant 
of time. Therefore, in order for a quasi-steady state to exist, the time derivatives in 
equations (3a) and (3b) must be negligible. This implies that the term (l 2 / vt ^(Pr/Ra) 1 / 2 
must be small. In order to evaluate this term, a characteristic time Tq, which is re- 
lated to the boundary conditions, must be chosen. Gebhart (ref. 7) shows that, for a ver- 
tical flat plate in an infinite medium, the characteristic time t q can be taken as 


2 

r 0 =— (bRa*) _2/5 (6) 

a 

where b is a function of Pr, and Ra* is the modified Rayleigh number. This same 
characteristic time may be associated with the problem of a confined fluid, since the wall 
temperature rise leads all other temperature rises in the system, and the rate of wall 
temperature rise is at least as great as the rate of rise of any other temperature in the 
system. The use of equation (6) produces the following condition under which a quasi- 
steady state may exist: 


(? rRa > 1/2 » 1 (7) 

(bRa*) 2 / 5 

Similarity Solution of Boundary Layer Equations 

It is assumed that a constant wall heat flux is imposed on a confined fluid in such a 
manner that a boundary layer flow exists and the system is in a quasi- steady state. If 
the interior temperature is of a certain form, the boundary layer temperature and veloc- 
ity profiles can be shown to exhibit the property of similarity. However, the interior 
temperature cannot be obtained directly from a solution of the interior equations, since 
no boundary conditions exist to govern these equations. Solutions for the boundary layer 
flow field can be generated if the interior equations are used to provide boundary condi- 
tions for the boundary layer equations. It has been shown in reference 5 that, in the in- 
terior of the fluid, 3T/3y « 0. If this relation is used in the interior equations, it follows 
that difz/dy « 0. Therefore, the following set of approximate boundary conditions can be 
written for the quasi- steady boundary layer equations: 
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( 8 ) 


The assumption that 3T/3y « 0 in the interior flow does not hold both near the liquid 
surface, where high values of the radial gradient have been observed, and near the tank 
bottom, where the boundary layer assumptions are invalid. 

It is possible to impose an artificial condition on the boundary layer equations by re- 
placing the term AT by the term T(x) - Tg(x) . (Note that the interior temperature T B 
is assumed to depend on x only.) This leads to the following partial differential equa- 
tions for the laminar boundary layer flow of a confined fluid in a quasi- steady state: 


'd}PJL-Mi JL.\ Mi = f/3(T - T r ) + v 
Jy dx dx dy) dy ay 3 


jl _ Ml JL\ t - a 

dy dx dx dy) 9y 2 


(9a) 


(9b) 


For the case of constant wall flux, the following similarity transformation may be per- 
formed on these equations. Let an independent variable be defined by 


n = 



( 10 ) 


where Gr* is the modified Grashof number. 

X 

New dependent velocity and temperature variables, F and 6, respectively, are de- 
fined by 


Ffa) = 


5 4 ' / 5 ^( Gr *) 1 / 5 


( 11 ) 
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Applying the similarity relations (10), (11), and (12) to equations (9a) and (9b) results 
in the following set of ordinary differential equations: 


F m + 4 FF m - 3F * 2 -0 = 0 

(13a) 

+ Pr[40'F - (Kj + 0)F*] = 0 

(13b) 


where the primes denote differentiation with respect to 77 . These equations hold only if 
the interior temperature is of the form 

T b (x) = T oo + AKjX 1 / 5 (14) 




Similarity variable, t? = (GRj/5) 1 ^ y/x 

Figure 2. - Dimensionless velocity and temperature 
distributions in boundary layer for flat plate in in- 
finite medium (ref. 8). 


where A and Kj are constants. For the 
present case of constant wall heat flux, A = 

The boundary conditions governing equa- 
tions (13a) and (13b) now become 


F = 0 
F' = 0 
0 * = 0 

F' = 0 

0 = 0 


77 = 0 


> 


>77 = 00 


(15) 


Equations (13a) and (13b) are coupled ordinary 
differential equations with Pr and Kj as 
parameters. Numerical solutions to these 
equations have been given in reference 8 for 
the case of = 0. The velocity and temper- 
ature profiles in the boundary layer for this 
case are shown in figure 2. These profiles 


8 




depend only on 77 for a given Pr; that is, the profiles exhibit the property of similarity. 
The case = 0 is the usual starting point for approximate propellant heating analyses. 


APPARATUS AND PROCEDURE 

A cutaway view of the test rig used for measuring boundary layer and interior tem- 
peratures is shown in figure 3. With the exception of the test section, the rig was the 
same as that used for the experiments discussed in reference 1. The test section (fig. 4) 
was a rectangular -parallelepiped glass tank with inside dimensions of 20. 3 centimeters 
wide by 30. 5 centimeters high by 5. 08 centimeters deep. The walls and bottom were 
fabricated from 0. 3 17 5- centimeter heat-resistant glass plate. A false bottom of phenolic 
resin (5. 08 cm by 20. 3 cm by 0. 635 cm) was located 3. 81 centimeters from the bottom 
of the test section and provided a definite starting place for sidewall heating. An outlet 
was located at the center of the bottom and false bottom plates. The top of the tank was 
open to permit the insertion of instrumentation. The sides of the test section from the 
false bottom to the top were painted black, while the bottom and the sides between the 
bottom and the false bottom were covered with heat-reflecting tape. The faces were 
clear to permit visual observations. 

Cylindrical infrared lamps, 30. 5 centimeters long by 0. 952 centimeter in diameter, 
were mounted in parabolic reflectors parallel to the blackened walls. The lamps pro- 
vided a wall heat flux that was proportional to the variable lamp power. 

The test fluid was a 2 to 1 mixture by volume of 1, 1, 1-trichloroethane and ethyl al- 
cohol. Table I summarizes the physical properties of these liquids and of the plate glass. 
Appendix B presents the method employed to obtain the thermal conductivity of the test 
liquid. An epoxy resin was used to cement the pieces of the test section together. 

Boundary layer temperature profiles were measured with two 0. 00 508- centimeter 
(2-mil) diameter copper -constantan thermocouples; a micrometer screw was used to 
position the temperature probe. The movable probe is shown in two positions in figure 3, 
one near the wall and the other on the tank centerline. Figure 4 shows a closeup view of 
the test tank and probe. Two automatic recording potentiometers were employed for all 
temperature measurements. 

Past experience had shown that end effects at the faces of the tank were negligible, 
as were the radial temperature gradients. Therefore, the tank acted as a two- 
dimensional test section, and only centerline and boundary layer temperatures were mea- 
sured. Centerline temperatures were measured with a rake which consisted of a 0. 952- 
centimeter-diameter by 30. 5- centimeter- long phenolic resin rod. Copper -constantan 
thermocouples were mounted in 1. 27-centimeter intervals starting at one end of the rod. 
The thermocouples extended 0. 3175 centimeter from the rod. 
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TABLE I. - PROPERTIES OF TEST MATERIALS 


Material 

Density, p, 

Specific heat, c p , 

Viscosity, p, 

Volume expansivity, /3, 

Thermal conductivity, k, 

Prandtl 

number, 

Pr, 

C pM- 

k > 
dimen- 
sionless 

kg 

lb mass 

J 

Btu 

(N) (sec) 

lb mass 

i 

°c 

i 

°F 

J 

Btu 

cu m 

cu ft 

(kg)(°C) 

(lb mass) (°F) 

sq m 

(hr) (ft) 

(sec)(m)(°C) 

(hr)(ft)(°F) 

Ethyl alcohol 

7.669X10 2 

0.478X10 2 

2.782X10 3 

0.665 

7.4X10' 4 

1.81 

1. 065xl0 -3 

0. 592X10" 3 

0.1547 

0.0895 

13.3 

(c 2 h 5 oh) 












1,1, 1-Trichloroethane 

1.300 

.812 

1.075 

.257 

6.0 

1.464 

1. 533 

.852 



— 

(CH 3 CC1 3 ) 












a Solution of 

1.122 

.701 

1.461 

.349 

6.4 

1.578 

1.377 

.765 

b . 0634 

.0367 

15.0 

two parts 1, 1, 1- 












trichloroethane 












and one part 












ethyl alcohol 





i 







Heat-resistant glass 





— 


J 


1.1311 

.654 

— 


a Ref. 9 was used to calculate properties of the mixture when applicable. 
b See appendix B. 




Figure 5. - Qualitative plot of temperature rise as function of time for wall and 
interior temperatures. Period during which system is in non-quasi- steady 
state, 0 to Tj? time interior temperature begins to rise, Tj; for x>Tj system 
is in quasi-steady state; average time during which first set of boundary layer 
temperatures is measured, (x 2 + t 3 )/2 = Taj; average time during which second 
set of boundary layer temperatures is measured, (x 4 +x 5 )/2 = x a2 . 

To measure temperature profiles near the wall, the temperature probe (see fig. 4) 
was positioned with one thermocouple in contact with the wall and the other approximately 
3 centimeters from the wall. The tank was filled to a height of 12.7 centimeters, and 
power was applied to the infrared lamps, which were then permitted to operate for ap- 
proximately 180 seconds before data were taken. Temperature readings were then taken 
simultaneously with both thermocouples, as the probe was moved in 0. 002 54- centimeter 
(1-mil) increments away from the wall. The temperature recorded by the thermocouple 
farthest from the wall was used as the temperature in the interior of the fluid Tg(x) . 

A qualitative plot of the wall and interior temperatures as functions of time is shown 
in figure 5. The wall temperature increases rapidly at first and the system is in a non- 
steady state. At some time t = t^, the interior temperature begins to rise, and the sys- 
tem approaches a quasi- steady state. Boundary layer temperature measurements were 
made at times greater than r^. Obviously, the measurements for a given value of x 
and q w extended over some period of time, for example, T 2 to r^. The average time 
after the start of heating t & during the measurement is given by (rg + Xg)/2. If a second 
measurement of the boundary layer temperature profile is made during a later time inter- 
val to Tg, the average time after the start of heating is given by (x^ + r g)/2 for this 
measurement. 

To measure the wall heat flux produced by the lamps, a second test tank was con- 
structed with thermocouples embedded on both sides of one wall. This test tank was 
identical to that used in reference 1 and is shown schematically in figure 6. This figure 
also shows the instrumented wall with the 16 copper -constantan thermocouples in place; 
the thermocouples were positioned opposite each other, 8 on either side of the wall. 

The inner thermocouples (those nearest the liquid) were mounted by drilling at an angle 
through the glass wall from the outside, inserting the thermocouple, and filling the re- 
maining space with an epoxy resin. The outer thermocouples were positioned on the 
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Figure 6. - Test section with instrumented wail for measuring wall heat flux. (All dimensions 
in centimeters. ) 


glass and held in place with epoxy. A small sheet of metal was placed over, but not in 
contact with, each of the outer thermocouples to prevent direct absorption of the radia- 
tion. The outside wall was then painted black. Temperature gradients through the wall 
were measured at six positions for two values of lamp power. The heat flux q w was 
then calculated from the equation 


"w 


= -k 


AT 
g AY 


where AY is the wall thickness (0. 3175 cm) and AT is the temperature difference 
across the wall. It was assumed that the wall temperature could be approximated as a 
linear function of Y. An experimentally determined relation between the lamp power P 
and the wall heat flux was found to be q w j/q w 2 = ^ 1 ^ 2 ’ w ^ ere the subscripts denote 
different values of applied lamp voltage. 

Figure 7 presents a plot of q at six positions on the wall for a value of P of 

w 

256 watts. Data are presented in 60-second intervals up to a maximum of 480 seconds 
after the start of heating. It can be seen that q w rises rapidly for 180 seconds and then 
remains fairly constant at any position. There were, however, differences between the 
values of q w measured at the various positions for a given time. These deviations were 
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TABLE H. - AVERAGE VALUES OF 



Figure 7. - Measured wall heat flux imposed by infrared lamps. 
Applied lamp power, 256 watts. 


WALL HEAT FLUX 
[Lamp power = 256 W. ] 


Time after start 
of heating, 
t , 
sec 

Average wall heat flux, q w , 

W/sq dm 

Btu/(hr)(sq in. ) 

60 

10. 89 

2. 40 

120 

13. 39 

2. 95 

180 

14.25 

3. 14 

240 

14. 48 

3. 19 

300 

13.62 

3.00 

360 

13. 89 

3. 06 

420 

14. 11 

3. 11 

480 

14.61 

3. 22 


caused by two experimental problems: first, the thermocouples were mounted with epoxy 
that could not be made uniformly thick; and, second, the lamps did not give a completely 
uniform distribution. 

The six values of q w were then arithmetically averaged for each time increment, 

and the average values are presented in table n. It can be seen that the average wall 

heat flux rises rapidly for 180 seconds and then remains fairly constant in tinie. Since 

all boundary layer temperatures were measured after at least 180 seconds of heating, the 

average of the average values of q_, 7 after 180 seconds was used as the value of q . For 

w w 

P equal to 256 watts, q w is 14. 16 watts per square decimeter (3. 12 Btu/(hr)(sq in.). 

The values of wall heat flux used in this experiment were chosen to give laminar 
boundary layer flow, which exists for values of Gr x up to 10 9 . Under this restriction, 
three values of q w were decided on, such that an arbitrary factor of 3 separated the 
maximum and minimum values. 


RESULTS AND DISCUSSION 


The underlying motive of these experiments was to determine whether boundary layer 
temperature profiles can be predicted by the similarity approach discussed in the analy- 
sis. The important parameter was found to be Kj, which describes the coupling between 
the interior and bulk temperatures. The region of interest was that part of the flow field 
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which was sufficiently far from the bottom and the liquid 
surface so that geometrical and free surface effects could 
be neglected. 



(a) Wall heat flux, 16.5 wafts per square 
decimeter (3. 69 Btu/(hr)(sq in. )). 



(b) Wall heat flux, 11.0 watts per square 
decimeter (2. 46 Btu/(hrMsq in. )). 



Temperature rise, °C 

iii. I.I.] 
0 10 20 30 40 

Temperature rise, °F 


(c) Wail heat flux, 5. 5 watts per square 
decimeter (1. 23 Btu/(hr)(sq in. )). 


Figure 8. - Centerline temperature pro- 
file with time after start of heating as 
parameter. Liquid level, 12. 7 centi- 
meters. 


Temperature Measurements 

Fluid temperatures taken on the central tank axis 
are presented in figure 8 for three values of q w> 5. 5, 

11. 0, and 16. 5 watts per square decimeter (1. 21, 2. 42, 
and 3. 63 Btu/(hr)(sq in.)), for times after the start of 
heating of 180, 300, 420, 540, and 660 seconds. Each of 
the profiles can be represented approximately by a func- 
tion of the form 


T B (x) = T co + K 2 x3/2 ( 16 ) 

Radial temperature profiles taken near a heated 
wall are presented in figure 9. Both wall heating and 
the boundary layer are assumed to start at the false bot- 
tom of the test section; therefore, the origin of the coor- 
dinate system is at the junction of the tank wall and the 
false bottom. Since measurements start at two times 
after the start of heating (see fig. 5, p. 13) for each 
value of x and since measurement of a single profile 
extends over a time interval of about 200 seconds, the 
average time after the start of heating r is used to 

a 

characterize each profile. Temperature distributions 
are given in terms of the following nondimensional vari- 
ables: 


T(x) - T b (x) 
* " V x > - T B< X > 
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Temperature variable, < p = (T-Tg)/(T W -Tg) 




(c) Average time after start of heating, 240 seconds; (PrRa) 1 ^ 2 = 
2. IxlO 4 , Grashof number, 2xl0 6 ; distance from start of 
boundary layer, 2. 54 centimeters. 

Figure 9. - Measured temperature distributions 




heated wall with wall heat flux as parameter. 
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Temperature variable, cp = (T-Tg)/(T W -Tg) 




(f) Average time after start of heating, 480 seconds; (PrRa) 1 ^ 2 = 
1. lxlO 5 ; Grashof number, 5. 8xl0 7 ; distance from start of 



(g) Average time after start of heating, 240 seconds, (PrRa) 1 ^ 2 = (h) Average time after start of heating, 480 seconds, (PrRa) 1/2 = 

2. 4xl0 5 ; Grashof number, 2. 5xl0 8 ; distance from start of 2. 3xl0 5 ; Grashof number, 2. 3xl0 8 - distance from start of 

boundary layer, 7. 62 centimeters. boundary layer, 7. 62 centimeters. 


Figured - Continued. 
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1. 21 

□ 
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0 
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Theoretical curve for constant 
interior temperature in liquid 
with Prandtl number of 15 


(i) Averaqe time after start of heating, 240 seconds; IPrRar' 
3.4xUn Grashof number, 5. 3x10* distance from start of 
boundary layer, 10. 16 centimeters. 


(j) Average time after start of heating, 480 seconds; (PrRar 
3. 3x10 ,• Grashof number, 4. 8xl0 8 ; distance from start of 
boundary layer, 10. 16 centimeters. 






.4 .8 1.2 1.6 0 

Similarity variable, 17 = (GrJ/5) 

(k) Average time after start of heating, 240 seconds; (PrRa) 1/2 = (I) Averag 

2. 5xl0 5 ; Grashof number, 2. 7xl0 8 ; distance from start of 2. 4xl0 5 ; 

boundary layer, 11. 43 centimeters. boundar 


(I) Average time after start of heating, 480 seconds; (PrRa) 1 

2. 4xl0 5 ; Grashof number, 2. 6xl0 8 ; distance from start of 
boundary layer, 11. 43 centimeters. 


Figured - Continued 




0 .4 .8 1.2 1.6 0 .4 .8 1.2 1.6 

Similarity variable, r\ » (Gr£/5) 1/5 y/x 


(m) Average time after start of heating, 240 seconds; (PrRa) 1/2 = (n) Average time after start of heating, 480 seconds; (PrRa) 1/2 * 

2. 6xl0 5 ; Grashof number, 3xl0 8 ; distance from start of 2. 53x10^- Grashof number, 2. 9xl0 8 ; distance from start of 

boundary layer, 12.065 centimeters. boundary layer, 12. 065 centimeters. 

Figured - Concluded. 

The temperature distributions with q w as a parameter are presented in figure 9. 

Distributions are given for two values of t , 240 and 480 seconds. As was shown in the 

a 

section GOVERNING PARAMETERS, boundary layer flow exists along the wall when the 

Rayleigh number is large or when (PrRa)"^ » 1. This quantity is given in each figure 

and indicates when boundary layer flow exists. The local Grashof number Gr x gives an 

indication of the flow regime within the boundary layer. The boundary layer flow is lami- 

9 

nar for Grashof numbers up to 10 , and, as can be seen from the plots, all the data fall 
in this regime. In each case, the theoretical profile obtained from reference 8 for con- 
stant interior temperature is shown along with the data. 

Temperature profiles for an x of 0. 238 centimeter (figs. 9(a) and (b)) show poor 
agreement with theory. This behavior is expected since the theoretical curve was de- 
rived by using the boundary layer assumptions and, therefore, is not valid near the start 
of the boundary layer. Figures 9(c) to (f) present the profiles for values of x of 2. 54 
and 5. 08 centimeters. Although there is considerable scatter in the data, experiment 
and theory appear to agree. The scatter is probably due to the small differences between 
the wall and bulk temperatures in this region and the inability of the experimental system 
to distinguish these differences accurately. The data for values of x of 7. 62, 10. 16, 
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and 11.43 centimeters (figs. 9(g) to (&)) 
agree very well with the theoretical pro- 
files for rj near zero. However, the 
experimental values of cp drop to zero 
at lower values of tj than do the theoret- 
ical values; this effect becomes greater 
with increasing x. The profiles for x 
of 12. 065 centimeters are shown in fig- 
ures 9(m) and (n). For this case, cp 
drops to zero very quickly and there is 
only poor agreement between theory and 
data even for values of 17 near zero. 


0 


.4 .8 

Boundary layer thickness, 2 

Figure 10. - Comparison of measured boundary layer shape with that 
predicted for liquid with constant interior temperature. 


Boundary Layer Characteristics 

For purposes of this report, the 
thickness of the temperature boundary 
layer will be defined as the value of r\ 
when <p has a value of 0. 2. Figure 10 
shows the shape of the boundary layer 

as determined from the data of figure 9. Also shown in the figure is a sketch of the 
boundary layer along a flat plate in an infinite medium (constant interior temperature). 

As x increases, the effects of stratification and the liquid surface become predominant, 
and the experimental boundary layer thickness approaches zero near the liquid surface. 

Heat is transferred from the walls to the interior of the fluid through the boundary 
layer. Therefore, the heat transfer coefficient between the wall and the boundary layer 
is useful in calculating interior temperature rises. As shown in the section GOVERNING 
PARAMETERS (see eq. (3c)), the local Nusselt number can be used to represent the local 
heat transfer coefficient. Figure 11 presents the experimental variation of Nu„ with the 
modified Grashof number 

infinite medium. The modified Grashof number Gr. 
layer length and is not indicative of the wall condition. The agreement between theory 
and data is quite good for values of x up to and including 10. 16 centimeters. This is 
reasonable since the slope of the boundary layer temperature profile evaluated at the wall 
determines Nu x , and it was indicated in the discussion of the temperature measurements 
that the experimental data and the theoretical curves agree quite well at r] = 0. For val- 
ues of x of 11. 43 and 12. 065 centimeters, Nu x is much larger than predicted because 
of the influence of the liquid surface. 


Figure 11 presents the experimental variation of Nu x 
Gr* along with the theoretical variation for a flat plate in an 

X 

* is used as a measure of boundary 
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TTHII I 

Wall heat flux, 

q w 

W/sq dm Btu/(hr)(sq in.) 


o 

5.5 

1.23 

□ 

11.0 

2.46 

o 

16.5 

3.69 


Theoretical curve for flat plate 

in infinite medium with 
Prandtl number of 15 

Plain symbols denote data measured at 
distances greater than 1. 27 cm from 
liquid surface. 

Tailed symbols denote data measured 
within 1. 27 cm of liquid surface 
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Figure 11. - Comparison of experimental data with predicted variation of Nusselt number with modified Grashof number. 


TABLE IH. - EXPERIMENTAL PARAMETERS 


FOR QUASI- STEADY OPERATION 


Wall 

heat 

flux, 

V 

W/sq dm 

Average time 
after start of 
heating, 

T a ? 

sec 

Average 
difference 
between wall 
and bulk tem- 
peratures, 

AT o’ 

°C 

Quasi- steady 
state term, 

(PrRa) 1 / 2 

(bRa*) 2/5 

5. 5 

240 

4. 64 

53. 3 

5. 5 

480 

4. 43 

52. 0 

11. 0 

240 

7. 16 

50. 0 

11.0 

480 

7. 16 

50. 0 

16. 5 

240 

9. 06 

48.0 

16. 5 

480 

8. 00 

45. 1 


Comparison of Data With Analysis 

The condition for quasi- steady opera- 
tion (eq. (7)) depends, through the Ray- 
leigh number, on the difference between 
the wall and interior temperatures at a 
given axial height. After 180 seconds of 
operation, this temperature difference 
was found to be nearly constant; that is, 
the wall and interior temperatures rise at 
approximately the same rate. Table HI 
shows the average temperature difference 
between the wall and interior, along with 
the value of the quasi- steady state term 
(PrRa) 1//2 (bRa*)' 2 / 5 , after 240 and 480 
seconds of operation for each value of 
heat flux employed. The maximum value 
of x is used to evaluate the quasi- steady 
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Temperature variable, <p = (T - Tg)/(T W - Tg) 



0 . 4 .8 1.20 .4 .8 1.2 

(a) Wall heat flux, 5. 5 watts per square deci- (b) Wall heat flux, 11. 0 watts per square deci- 
meter. meter 



0 .4 .8 1.2 


Similarity variable, r? = (Gr£/5)^ y/x 

(c) Wall heat flux, 16. 5 watts per square deci- 
meter. 

Figure 13. - Measured temperature distribution near heated wall with distance from start of boundary 
layer as parameter. Average time after start of heating, 240 seconds. 
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condition, since this condition depends on The results show that a quasi- steady 

state existed for values of r of 240 and 480 seconds. 

The temperature profiles in the boundary layer possess the property of similarity if 
they depend only on the similarity variable rj for a given heat flux. Figure 12 shows 
measured boundary layer temperature profiles for a q w of 5. 5 watts per square deci- 
meter with r a as a parameter for values of x of 5. 08, 7. 62, 10. 16, and 11. 43 centi- 
meters. The data are independent of t & for these values of x. The same behavior was 
noted for the values of q w not shown in figure 12. 

Profiles with x as a parameter are presented in figure 13. Data are presented for a 
q w of 5. 5, 11. 0, and 16. 5 watts per square decimeter for a t & of 240 seconds. The 
data indicate that the profiles are approximately independent of x for x-values of 5. 08, 
7.62, and 10. 16 centimeters. This is true, even though Tg(x) is not of the required 
form. The nondependence of the profiles on t_ and x between 5. 08 and 10. 43 centi- 

a 

meters indicates that the concept of similarity can be used in this region. 

For values of x greater than 10. 43 centimeters, similarity does not exist. The ap- 
parent cause is that the liquid surface perturbs the boundary layer flow and destroys 
similarity characteristics in this region. This perturbation is reflected in the local 
Grashof numbers, which decrease for values of x greater than 10. 43 centimeters (see 
figs. 9(k) to (n)); this effect indicates that the difference between the wall and interior 
temperatures is decreasing as the surface is approached. 

Equation (13b) shows that the interior and boundary layer temperatures are coupled 
through the parameter Kj. This equation predicts that the boundary layer temperature 
profiles will possess the characteristic of similarity provided that the interior tempera- 
ture is of the form shown in equation (14). Equation (13b) also shows that the boundary 
layer temperature profiles will be nearly identical to those for constant interior temper- 
ature if Kj « [0 (tj) |. It has been found experimentally that the interior temperature is 
not of the required form (see eq. (16)). Nevertheless, the boundary layer temperature 
profiles appear to be similar in certain regions of the boundary layer, and the profiles 
for Kj = 0 give a good approximation of the data. It is of interest to express these ideas 
in a quantitative manner. 

Since the experimental boundary layer temperature profiles exhibit similarity to 
some degree, it can be assumed that the interior temperature profile is of the required 
form. This means that 


AKjX 1 / 5 = K 2 x 3 / 2 (17) 

Because A and Kg are assumed to be constants, this equation infers that Kj is actu- 
ally a function of x. To find an explicit form for K^(x/L), divide both sides of equa- 
tion (17) by L and solve for Kj. Then, 
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*1 = K 2 L 


(18) 


13/10 ^ 13 / 10 A -1 


Since Kj is a function of x, the boundary layer profiles cannot be truly similar, and it 
was noted earlier in this section that the profiles were only approximately independent 
of x. The slight x dependence was most obvious for 77 » 0. Now 0 ( 77 ) has a maximum 
absolute value when 77 = 0 and approaches a value of zero as 77 — <*>, that is, in the inte- 
rior of the fluid. Assume that Kj « 1 0 (0) | . Then, as rj increases, a point is reached 
at which = 1 0 (77 ^) | . In the region near 77 the profiles cannot be similar. The posi- 
tive ratio Kj/[-0(r7)] can be formed by using equations (12) and (17). Then, 


K^x) ^ t 
-0(77) ~ AT 1 
AT B 


(19) 


Therefore, Kj(x) = 1 0 (77 ^) | when AT = 2 ATg. In plotting the data, the parameter <p(i 7 ) 
given by <p = 0(r 7 ) /0 (0) was used. Then cpfyi^) can be found by setting AT = 2 ATg and 


<p(Vl) = 2 ATg 


AT B 

AT w " at b 


or 


<pfa 1 ) = ( 20 ) 

AT W t 

aT B 

This equation is somewhat ambiguous because AT W and ATg change with 77 . How- 
ever, values of these quantities when, for example, 77 = 1 . 0 can be used to provide a 
close approximation for cpfy}^)- For the system under consideration, it was found that 
AT^y » ATg when 77 = 1.0, and from equation (20), (pfy]^) is nearly zero. This agrees 
with the observation that the profiles are similar when cp is near 1 . 0 , but similarity is 
less evident as <p approaches zero. 

To make the foregoing discussion more concrete, consider the case where q w = 16. 5 
watts per square decimeter and r = 660 seconds (fig. 8 , p. 16). For this case, 

ATg = 0. 40 x 3 / 2 . Then K 2 = 0. 40, A = 88 . 79, and l 13//1 ° = 0. 0685 in the mks system 
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of units. A few values of Kj(x/L) as calculated from equation (18) would be x/L = 0. 2 
and Kj = 6xl0" 5 , x/L = 0. 5 and Kj = 1. 5xl0" 4 , and x/L = 1. 0 and Kj = 3xl0" 4 . Thus, 
Kj is of the order of 10” 4 for the case being considered. When equation (12) and the 
measured temperature difference Tg - T^ are used, 0(0) is of the order of 5x10" 4 for 
those values of x/L where similarity was observed. Furthermore, 77 ^ is very close to 
1.0 for this case, and 8 ( t ] j ) is of the order of 10~ 4 . 

CONCLUDING REMARKS 

The analysis and experimental results show that under the conditions studied, that is, 
for constant wall heat flux on a confined liquid of Prandtl number 15, the following are 
true: 

1. The flow fields are in a quasi-steady state when the relation (PrRa) 4 ^^(bRa*)"^^^ 
» 1 is satisfied. 

2. There is a region in the laminar boundary layer in which the temperature profiles 
exhibit the property of similarity and are predicted fairly well by the infinite medium, 
flat plate model. 

3. There is a region in the boundary layer near the liquid surface in which surface 
effects destroy similarity characteristics and reduce the boundary layer thickness. 

4. The coupling between the boundary layer and interior temperature profiles is 
weak if [k 1 /|0(O)|] « 1. 

(The foregoing conclusions are based on the assumption that boundary layer flow exists 
when (PrRa) 4/ ^ » 1.) 

These conclusions, however, are based on temperature profiles only. To complete 
the study of the laminar boundary layer in a confined liquid, velocity measurements 

Q 

should also be obtained. Furthermore, the transition and turbulent regimes (Gr x > 10 ) 
must be investigated in order to obtain a complete understanding of the boundary layer in 
an enclosed liquid. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 19, 1966. 


27 



APPENDIX A 


SYMBOLS 


A 

constant, 

Pr 

Prandtl number, v/a, dimension- 

B 

body force per unit volume, 


less 


N/cu m 

P 

pressure, N/sq m 

b 

constant, 0. 0090 for liquid with 

Po 

reference pressure, N/sq m 


Prandtl number of 15 

q w 

wall heat flux, W/sq dm 

C P 

specific heat, J/(kg) (°C) 

Ra 

Rayleigh number, Gr x Pr, dimen- 

D 

viscous dissipation function 


sionless 

F 

similarity function for velocity 

Ra* 

modified Rayleigh number, Gr*Pr, 

f 

body force, N 


dimensionless 

Gr x 

Grashof number based on x, 

T 

temperature, °C 


(x^f/3 AT)/u^, dimensionless 

t b 

bulk temperature, function of x, 

Gr* 

modified Grashof number based 


°C 

A 

on x, (fpq^x^/kv^, dimen- 
sionless 

Tco 

constant reference temperature, 
°C 

h 

heat transfer coefficient, 
J/(sec) (m) 2 (°C) 

AT 

T - T b , °C 

K 1 

function describing coupling be- 

AT 0 

reference temperature, °C 


tween boundary layer and inte- 
rior flows 

u 

x-component of boundary layer ve- 
locity, cm/sec 

K 2 

constant 

V 

y- component of boundary layer ve- 

k 

thermal conductivity of test 


locity, cm/sec 


fluid, J/(sec)(m)(°C) 

X 

coordinate parallel to heated wall, 

k 

s 

thermal conductivity of heat re- 


cm 

o 

sistant glass, J/(sec)(m)(°C) 

AY 

wall thickness, cm 

L 

liquid level, used as character- 
istic length, cm 

y 

coordinate perpendicular to heated 
wall, cm 

Nu x 

Nusselt number, (hx)/k, dimen- 
sionless 

a 

thermal diffusivity, k/pc^, 
sq m/sec 

P 

lamp power, W 

0 

volume expansivity, 1/°C 
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V similarity variable 

9 similarity function of temperature 

dimensionless temperature vari- 
able, AT/ATq 

M viscosity, N(sec)/sq m 

v kinematic viscosity, p/p, 

sq m/sec 

p density, kg/cu m 

t time, sec 


r average time after start of heat- 

a 

ing, sec 

Tq characteristic time, sec 

cp temperature variable, [T(x) - 

Tb( x )]/[ t w ( x ) - T b (x)], dimen 
sionless 

ip stream function, sq m/ sec 

Subscripts: 
i interior 

w wall 
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APPENDIX B 


CALCULATION OF THERMAL CONDUCTIVITY OF TEST FLUID 


The following definitions are made: 

thermal conductivity of ethyl alcohol 
kg thermal conductivity of trichloroethane 

kgj thermal conductivity of test fluid (2 to 1 mixture by volume of 1, 1, 1-trichloroethane 
and ethyl alcohol) 

k thermal conductivity of heat resistant plate glass 
& 

Although a fair amount of data exists on the values of kj and k g , no data could be 
found for the value of kg. Furthermore, there is no universally acceptable way of pre- 
dicting the thermal conductivity of an organic liquid or of a mixture of organic liquids. 
Consequently, the following scheme was devised for the calculation of kgj. 

The heat flux at the heated wall is given by the relations 


<L 


w 



9T 



% = -k 21 


5T 

3y 


w 


21 


where the subscripts indicate that the temperature gradient must be measured in the 
glass or in the liquid. It will be noted that q w was calculated in the section APPARATUS 


AND PROCEDURE by using the first of the aforementioned two relations 
9T 

9y 


Furthermore, 

was measured directly. With these data, kg j was obtained from the relation 


w 


21 


k 21 " “ 


a w 


9T 

3y 


w 


21 


Since k 21 is a function of T, values of k 21 were calculated for each value of q w 
and r & used for the three values of x of 5. 08, 7.62, and 10. 43 centimeters. All values 
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obtained in this manner were then arithmetically averaged to give kgj for this set of 
experiments. Table IV gives the calculated in this manner. 


TABLE IV. - CALCULATED VALUES OF 
THERMAL CONDUCTIVITY 


[Average k 21 = 0. 0634 J/(sec)(m)(°C) = 0. 0367 Btu/(hr)(ft)(°F). ] 


Wall 

Time after 

Distance 

Thermal conductivity, k 21 , 

heat 

start of 

from start 



flux, 

heating, 

of boundary 

J/(sec)(m) (°C) 

Btu/(hr)(ft)(°F) 

V 

sec 

layer, 



W/sq dm 


x, 





cm 



5. 5 

240 

5. 08 

0. 0667 

0. 0386 



7.62 

. 07 59 

. 0439 



10. 16 

. 0764 

. 0442 

5. 5 

480 

5. 08 

0. 0574 

0. 0332 



7. 62 

. 0622 

. 0360 



10. 16 

. 0511 

. 0296 

11.0 

240 

5. 08 

0. 0669 

0. 0387 



7.62 

. 0657 

.0380 



10. 16 

. 0537 

.0311 

11.0 

480 

7.62 

0. 0584 

0.0338 



10. 16 

. 0586 

. 0339 

16. 5 

240 

7. 62 

0. 0551 

0. 0319 



10. 16 

. 0733 

. 0424 

16. 5 

480 

5. 08 

0. 0655 

0. 0379 



7. 62 

. 0624 

. 0361 



10. 16 

. 0665 

. 0385 
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